Precision measurement of the 7 Be solar neutrino interaction rate in Borexino 
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The rate of neutrino-electron elastic scattering interactions from 862 keV 7 Be solar neutri- 
nos in Borexino is determined to be 46.0±1.5(stat)j^'g(syst) counts/(day-100ton). This corre- 
sponds to a ^-equivalent 7 Be solar neutrino flux of (3.10T0.15) x 10 9 cm _2 s _1 and, under the 
assumption of v a transition to other active neutrino flavours, yields an electron neutrino sur- 
vival probability of 0.51±0.07 at 862 keV. The no flavor change hypothesis is ruled out at 5.0 a. 
A global solar neutrino analysis with free fluxes determines <£ > p P =6.06to'Q 2 i x 10 10 cm~ 2 s _1 and 
< l > Civo<1.3x 10 9 cm _2 s _1 (95% C.L.). These results significantly improve the precision with which 
the MSW-LMA neutrino oscillation model is experimentally tested at low energy. 

PACS numbers: 13.35.Hb, 14.60. St, 26.65.+t, 95.55.Vj, 29.40.Mc 



In the past 40 years, solar neutrino experiments [IHl] 
have revealed important information about the Sun [5] 
and have shown that solar neutrinos undergo flavour 
transitions that are well described by Mikheyev-Smirnov- 
Wolfenstein Large Mixing Angle ("MSW-LMA") type 
flavour oscillations [§]. Reactor antineutrino measure- 
ments [7] also support this model. The MSW model 
predicts a transition in the solar v e survival probabil- 
ity ("P ee ") at neutrino energies of about 1-4 McV. This 
transition is currently poorly tested. Therefore, in order 
to test MSW-LMA more thoroughly, to probe other pro- 
posed neutrino oscillation scenarios [8j, and to further 



improve our understanding of the Sun, it is important 
that experimental measurements of the low energy solar 
neutrino fluxes be improved [Sj. At 862 keV, the abun- 
dant, mono-energetic, 7 Be solar neutrinos can provide a 
precise probe of the survival probability in this interest- 
ing region. In addition, a precise determination of the 
7 Be flux combined with existing results from radiochem- 
ical experiments [TJ [2] yields improved constraints on the 
pp and CNO solar neutrino fluxes. 

The Borexino experiment at Gran Sasso [TO] detects 
neutrinos through the neutrino-electron elastic scattering 
interaction on a ^278 metric ton liquid scintillator target. 
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The low energy backgrounds in the detector have been 
suppressed to unprecedented levels [11] , making Borexino 
the first experiment capable of making spectrally resolved 
measurements of solar neutrinos at energies below 1 MeV. 
We have previously reported a direct measurement of the 
7 Be solar neutrino flux with combined statistical and sys- 
tematic errors of 10% [T^]. Following a campaign of de- 
tector calibrations and a 4-fold increase in solar neutrino 
exposure, we present here a new 7 Be neutrino flux mea- 
surement with a total uncertainty less than 5%. For the 
first time, the experimental uncertainty is smaller than 
the uncertainty in the Standard Solar Model ("SSM") 
prediction of the 7 Be neutrino flux [13j^| 

The new result is based on the analysis of 740.7 live 
days (after cuts) of data which were recorded in the pe- 
riod from May 16, 2007 to May 8, 2010, and which cor- 
respond to a 153.6 ton-yr fiducial exposure. 

The experimental signature of 7 Be neutrino inter- 
actions in Borexino is a Compton-likc shoulder at 
^660 keV. Fits to the spectrum of observed event ener- 
gies are used to distinguish between this neutrino scatter- 
ing feature and backgrounds from radioactive decays |12j . 
Two independent fit methods were used, one which is 
Monte Carlo based and one which uses an analytic de- 
scription of the detector response. In both methods, 
the weights for the 7 Be neutrino signal and the main 
radioactive background components ( 85 Kr, 210 Po, 210 Bi, 
and n C) were left as free parameters in the fit, while 
the contributions of the pp, pep, CNO, and 8 B solar neu- 
trinos were fixed to the SSM-predicted rates assuming 
MSW neutrino oscillations with tan 2 6» 12 =0.47±°;o4 and 
Am 2 2 =(7.6±0.2)xl0~ 5 eV 2 [H]. Th e impact of fixing 
these fluxes was evaluated and included as a systematic 
uncertainty. The rates of 222 Rn, 218 Po, and 214 Pb sur- 
viving the cuts were fixed using the measured rate of 
214 Bi- 214 Po delayed coincidence events. The Monte Carlo 
method also includes external 7-ray background, which 
makes it possible to extend the fit range in this method 
to higher energies. The energy scale and resolution were 
floated in the analytic fits, while the Monte Carlo ap- 
proach automatically incorporates the simulated energy 
response of the detector. 

The stability of each fit method was studied by repeat- 
ing the fits with slightly varied fit characteristics (e.g. 
fit range and histogram binning) and different methods 
of data preparation. The latter included changing the 
method used to estimate the event energies, and varying 
the pulse shape analysis ("PSA") technique [TS] used to 
remove 210 Po and other a events between a highly effi- 



1 Throughout this Letter we use the high metallicity SSM predic- 
tions from the "GS98" column in Table 2 of |13| as our reference 
SSM. For comparison, the 7 Be neutrino flux predicted by the low 
metallicity model (also given in [13]) is 8.8% lower than the high 
metallicity prediction. 



Fit: %7NDF = 141/138 
7 Be: 45.5 ± 1.5 
85 Kr: 34.8 ± 1.7 
210 Bi: 41.5 ± 1.5 
"C: 28.9 ±0.2 
210 Po: 656.0±9.8 
External: 4.5 ± 0.7 
pp, pep, CNO (Fixed) 
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Fit: Z 2 /NDF = 99/95 




7 Be: 47.0 ± 1.9 




85 Kr: 24.6 ± 3.2 




— 210 Bi: 40.6 ±2.6 
"C: 28.0 ±0.4 




pp, pep, CNO (Fixed) 
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FIG. 1. Two example fitted spectra; the fit results in the leg- 
ends have units [counts/(day-100ton)]. Top: A Monte Carlo 
based fit over the energy region 270-1600 keV to a spectrum 
from which some, but not all, of the a events have been re- 
moved using a PSA cut, and in which the event energies were 
estimated using the number of photons detected by the PMT 
array. Bottom: An analytic fit over the 290-1270 keV energy 
region to a spectrum obtained with statistical a subtraction 
and in which the event energies were estimated using the to- 
tal charge collected by the PMT array. In all cases the fitted 
event rates refer to the total rate of each species, independent 
of the fit energy window. 

cient statistical subtraction method [H] and a cut-based 
technique which removes a fraction of the a events with 
a very small loss of /3 events. The example spectra shown 
in Fig. [T] illustrate the stability of our fit procedure; the 
8 B neutrino and 214 Pb, 222 Rn, and 218 Po background 
spectra are small on the scale of the plots and are not 
shown. The results of these and other fits using different 
permutations of the fit characteristics and data prepara- 
tion techniques described above were averaged to obtain 
the central values reported in Table|l] the spread between 
the results is included in the systematic uncertainty. 
TABLE I. Average Fit Results [counts/(day-100ton)]. 



'Be 46.0±1.5(stat)+^(syst) 



85 Kr 31.2±1.7(stat)±4.7(syst) 
210 Bi 41.0±1.5(stat)±2.3(syst) 
n C 28.5±0.2(stat)±0.7(syst) 
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TABLE II. 7 Be Systematic Uncertainties [%]. 

Source [%\~ 
Trigger efficiency and stability <0.1 

Live time 0.04 

Scintillator density 0.05 

Sacrifice of cuts 0.1 

Fiducial volume —1'.3 

Fit methods 2.0 

Energy response 2.7 

Total Systematic Error —s'e 



The main systematic uncertainties in our measurement 
of the 7 Be interaction rate are listed in Table [n] The 
dominant contributions come from the determination of 
the fiducial volume, our understanding of the detector 
energy response, and the variation between the results 
of the different fit procedures. We note the significant 
decrease in the uncertainties associated with the detector 
energy response and the definition of the fiducial volume, 
by factors of 4.6 and 2.2 respectively, relative to [12]. This 
improvement was made possible by several campaigns of 
detector calibration. 

During the first 1.6 years of Borexino operation, the en- 
ergy and position reconstruction algorithms were tuned, 
and their performances estimated, using intrinsic activi- 
ties such as 14 C, 210 Po, and 11 C. The first deployed source 
calibrations were carried out in 4 campaigns between 
October 2008 and July 2009. Encapsulated radioactive 
sources, including 57 Co, 139 Ce, 203 Hg, 85 Sr, 54 Mn, 65 Zn, 
40 K, 60 Co, and 222 Rn, were placed inside the scintillator 
volume using a rod-based source deployment system. Us- 
ing seven CCD cameras mounted within the PMT array, 
the positions of the sources could be determined with a 
precision better than 2 cm. In aggregate, more than 35 
live-days of calibration data were recorded with sources 
at more than 250 positions within the scintillator volume. 

The systematic uncertainty in the definition of the fidu- 
cial mass was evaluated by comparing the position of the 
source as measured by the CCD camera system and as 
reconstructed using the PMT array for a number of runs 
with sources deployed near the boundary of the fiducial 
volume. The 7-ray sources were also used to validate 
and improve both the Monte Carlo (the simulated and 
observed optical responses agree at the 1.5% level within 
the fiducial volume) and the detector energy response 
function used in the analytic fitting procedure. The un- 
certainties in these tunings and in our understanding of 
the calibration data were included in the energy scale 
systematic. A dominant contribution to the latter came 
from the uncertainty in the scintillator quenching model 
used to extrapolate the detector optical response from 
the 7-ray events used for energy calibration to the single 
electron events which comprise the majority of the data. 

It may be noted in Table [i] that the 85 Kr rate has a 
larger systematic uncertainty than the other fitted rates. 
This is due to a larger variation in the 85 Kr rate between 



the different fit procedures. We typically obtain larger 
85 Kr values when fewer a events are removed from the 
spectrum before fitting and when the fit range is extended 
to encompass lower energy events. However, as the Be 
rate is mostly constrained by the spectral region between 
550 and 750 keV, the variation in the 85 Kr rate reflects 
only weakly in the 7 Be rate. We note that the fitted 
85 Kr rate is consistent with an independent measurement 
of 30.4±5.3(stat)±1.3(syst) counts/ (day-100ton) for the 
85 Kr activity obtained using 85 Kr- 85m Rb delayed coinci- 
dences. 

Our best value for the interaction rate of 
862 keV 7 Be solar neutrinos in Borexino is 
46.0±1.5(stat)t^(syst)counts/(day-100ton). If the 
neutrinos are assumed to be purely v e , this corre- 
spond^] to an 862 keV 7 Be solar neutrino flusj^] of 
(2.78±0.13)xl0 9 cm" 2 s" 1 . The corresponding flux 
prediction from the SSM is (4.48±0.31) x 10 9 cm^s" 1 , 
which, if all the neutrinos remained v ei would yield 
an interaction rate of 74.0±5.2 counts/(dayT00 ton) in 
Borexino; the observed interaction rate is 5.0 a lower. 
The ratio of the measured to the predicted ^-equivalent 
flux is 0.62±0.05. Under the assumption that the reduc- 
tion in the apparent flux is the result of v e oscillation 
to or v T (which undergo electron elastic scattering 
interactions, but with a cross section about 4.5 times 
lower than v e at this energy), we find P ee =0.51±0.07 
at 862 keV. The improved constraint on the low energy 
solar P ee is shown in Fig. [2] 

We have also performed a global analysis to determine 
the MSW neutrino mixing parameters in the two-flavour 
approximation. This included the rate and spectrum in- 
formation from the other solar neutrino experiments [T|- 
H] , the SSM flux predictions [13] (with the exception that 
the 8 B flux was left free), and the current result. We find 
best fit oscillation parameters of tan 2 6* 12 =0.468^q 030 
and Am 2 2 = (5.2+ ( l ) ;g)xl0" 5 eV 2 ; including the Kam- 
LAND reactor anti-neutrino data [7] these become 
tan 2 6»i 2 =0.457t 9 °g and Am 2 2 =(7.50t°;^) x 10" 5 eV 2 . 
The new result slightly improves the precision with which 
the MSW mixing parameters can be determined. 

Alternatively, by assuming MSW-LMA solar neutrino 
oscillations, the Borexino results can be used to measure 
the 7 Be solar neutrino flux. Using the oscillation pa- 
rameters from P3], the Borexino result corresponds to a 
total 7 Be neutrino flux $7 Bc =(4.84±0.24) x 10 9 cm" 2 s" 1 
(after oscillation effects, the SSM prediction corre- 
sponds to an 862 keV 7 Be neutrino interaction rate of 



2 In converting from interaction rates to fluxes, we use the electron 
scattering cross section from 1161 . with updated radiative correc- 
tion parameters from 1 141 1171 , and a scintillator electron density 
of (3.307±0.003)xl0 29 /ton. 

s Note the distinction between the " 7 Be solar neutrino flux" and 
the "862 keV 7 Be solar neutrino flux": the latter is an 89.6% 
branch of the former. 
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pp - all solar 
8 B - all solar 

All solar without Borexino 
MSW Prediction 
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FIG. 2. The global experimental constraints on the low energy 
solar P ee . For the 7 Be point, which shows the current result, 
the inner (red) error bars show the experimental uncertainty, 
while the outer (blue) error bars show the total (experimen- 
tal + SSM) uncertainty. The remaining points were obtained 
following the procedure in [Tj5], wherein the survival probabil- 
ities of the low energy (pp), medium energy, and high energy 
( 8 B) solar neutrinos are obtained, with minimal model de- 
pendence, from a combined analysis of the results of all solar 
neutrino experiments. To illustrate Borexino's effect on the 
low energy P ee measurements, the green (dashed) points are 
calculated without using the Borexino data. The MSW-LMA 
prediction is also shown for comparison; the band defines the 
1-ct range of the mixing parameter estimate in [14], which 
does not include the current result. 

47.5±3.4counts/(day-100ton) in Borexino). The ra- 
tio of our measurement to the SSM prediction gives 
/7 Bc =0.97±0.09. 

Finally, we have used our new result to update the 
global experimental constraint on the other low energy 
solar neutrino fluxes by performing a global analysis in 
which the 8 B, 7 Be, CNO, and pp fluxes were left as 
free parameters, following [19J . Under the luminosity 

(6.06±%?£)xl0 10 axT^- 1 and 

s 

fraction of the SSM predicted fluxes these correspond 
to / pp =1.013±g;^ and / CNO <2.5 at 95% C.L. The lat- 
ter limits the CNO contribution to the solar luminosity 
to <1.7% (95% C.L.). Both the precision of the pp flux 
determination and the constraint on the CNO flux are 
improved by approximately a factor of two by our new 
result. 

We have presented a measurement of the interaction 
rate of 862 keV 7 Be solar neutrinos in Borexino with a 
total uncertainty less than 5%. This precise measure- 
ment has also improved our experimental understanding 
of the other low energy solar neutrinos. The Borexino 
measurements have allowed us to test and validate the 
MSW-LMA model in the vacuum oscillation regime and 
at the lower energy edge of the transition region with 
unprecedented precision. 

The Borexino program is made possible by funding 
from INFN (Italy), NSF (USA), BMBF, DFG, and MPG 



constraint, we find & pp — ^.uu_ 06 ; 
$cno<1-3x10 9 cm _2 s _1 at 95% C.L. Expressed as a 
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